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Abstract 

We study phase transitions in models of opinion formation which are based on the 
social impact theory. Two different models are discussed: (i) a cellular-automata 
based model of a finite group with a strong leader where persons can change their 
opinions but not their spatial positions, and (ii) a model with persons treated as 
active Brownian particles interacting via a communication field. In the first model, 
two stable phases are possible: a cluster around the leader, and a state of social 
unification. The transition into the second state occurs for a large leader strength 
and/or for a high level of social noise. In the second model, we find three stable 
phases, which correspond either to a "paramagnetic" phase (for high noise and 
strong diffusion), a "ferromagnetic" phase (for small nose and weak diffusion), or a 
phase with spatially separated "domains" (for intermediate conditions). 
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1 Introduction 



During the last years there has been a great interest in applications of sta- 
tistical physics in social science [1,2]. Usually economical models are studied 
using the techniques of stochastic dynamics [3], percolation theory [4] or the 
chaos paradigm [5]. Another important subject of this kind is the process of 
opinion formation treated as a collective phenomenon. On the "macroscopic" 
level it can be described using the master equation or Boltzmann-like equa- 
tions for global variables [6-8], but microscopic models are constructed and 
investigated as well [9,10] using standard methods of statistical physics. A 
quantitative approach to the dynamics of opinion formation is related to the 
concept of social impact [11-17], which enables to apply the methods similar 
to the cellular automata approach [9,18]. 

The aim of this work is to study various kinds of phase transitions in two mod- 
els based on the social impact theory. In Sec. 2 we consider phase transitions 
in a social impact model that can occur in a finite group in the presence of a 
strong individual (a leader) [19-21]. As two special cases, we discuss a purely 
deterministic limit and a noisy model. Sec. 3 is devoted to an extension of 
social impact models to include phenomena of migration, memory effects and 
a finite velocity of information exchange. Here the concepts of active Brownian 
particles [22-24] and the communication field [25] will be applied. 



2 Phase transitions in the presence of a strong leader 

2.1 The model 

Our system consists of N individuals (members of a social group); we assume 
that each of them can share one of two opposite opinions on a certain subject, 
denoted as Oi = ±1, i = 1,2,...N. Individuals can influence each other, and 
each of them is characterised by the parameter > which describes the 
strength of his/her influence. Every pair of individuals (i,j) is ascribed a 
distance in a social space. The changes of opinion are determined by the 
social impact exerted on every individual: 

N 

where g{x) is an increasing function of social distance. (3 is a so-called self- 
support parameter reflecting the inclination of an individual to maintain his/her 
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current opinion, h is an additional (external) influence which may be regarded 
as a global preference towards one of the opinions stimulated by mass-media, 
government policy, etc. 



Opinions of individuals may change simultaneously (synchronous dynamics) 
in discrete time steps according to the rule: 



<Ji(t + 1) 



<Ji{t) with probability 
-di{t) with probability 



exp(-yr) 

exp(-VT) + exp(VT) 

exp(Ij/T) 

exp(-VT) + exp(/ l /T) 



(2) 



Eq. (2) is analogous to the Glauber dynamics with — ij <jj corresponding to the 
local field. The parameter T may be interpreted as a "social temperature" de- 
scribing a degree of randomness in the behaviour of individuals, but also their 
average volatility (cf [10]). The impact Ii is a "deterministic" force inclining 
the individual i to change his/her opinion if > 0, or to keep it otherwise. 
The model is a particular case of the system considered in [14]. 

We assume that our social space is a 2D disc of radius R >> 1, with the 
individuals located on the nodes of a quadratic grid. The distance between 
nearest neighbours equals 1, while the geometric distance models the social 
immediacy. The strength parameters Sj of the individuals are positive random 
numbers with probability distribution q(s) and the mean value s. In the centre 
of the disc there is a strong individual (who we will call the "leader"); his/her 
strength sl is much larger than that of all the others (sl 3> Sj). 



2.2 Deterministic limit 



Let us first recall the properties of the system without noise, i.e. at T = 
[19,20]. Then, the dynamical rule (2) becomes strictly deterministic: 

Oi(t + 1) = -sign(7j(7i). (3) 



Considering the possible stationary states we find either the trivial unification 
(with equal opinion ±1 for each individual) or, due to the symmetry, a circular 
cluster of individuals who share the opinion of the leader. This cluster is 
surrounded by a ring of their opponents (the majority). These states remain 
stationary also for a small self-support parameter (3; for sufficiently large (3 
any configuration may remain "frozen" . 

Using the approximation of continuous distribution of individuals (i.e. replac- 
ing the sum in (1) by an integral) one can calculate the size of the cluster, 
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i.e. its radius function of the other parameters. In the case of g(r) = r 

and s = 1 we get from the limiting condition for the stationarity / = at the 
border of the cluster: 



1 r 



2nR - VTri/3 -h±\l (2nR - ^ ± (3 - h) 2 - 32s L . (4) 




This is an approximate solution valid for a <C R, but it captures all the 
qualitative features of the exact one which can be obtained by solving a tran- 
scendent equation (cf. Fig. 1). Here and in the next section we assume that 
the leader's opinion is aj_, — +1, but the analysis is also valid for the opposite 
case if h — > —h. 

The branch with the "— " sign in front of the square root in Eq. (4) corresponds 
to the stable cluster. The one with "+" corresponds to the unstable solution 
which separates the basins of attraction of the stable cluster and unification 
(cf. Fig. 1). Owing to the two possible signs at the self-support parameter j3 
in (4), the stable and unstable solutions are split and form in fact two bands. 
The states within the bands are "frozen" due to the self-support which may 
be regarded as an analogy of the dry friction in mechanical systems. This way 
also the unstable clusters can be observed for j3 > and appropriately chosen 
initial conditions. 

According to Eq. (4) real solutions corresponding to clusters exist provided 



Otherwise the general acceptance of the leader's opinion (unification) is the 
only stable state. When, having a stable cluster, the condition (5) is violated 
by changing a parameter e.g. sl or h, one can observe a discontinuous phase 
transition: cluster — > unification. 

If, on the other hand, the leader's strength is too weak, it may be impossible for 
him/her not only to form a cluster but also to maintain his/her own opinion. 
The limiting condition for the minimal leader's strength simin to resist against 
the persuasive impact of the majority can be calculated from the limiting 
condition I L = (II - the impact exerted on the leader): 



Fig. 1 shows a phase-diagram s^-a for h = 0. All the plots are made for a space 
of radius R = 20 (1257 individuals) and /3 — 1 unless stated otherwise. Points 
in Fig. 1 are obtained by numerical simulations of (3) while the curves are 
solutions of a transcendent equation following from the stationary condition 




(5) 
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1(a) = 0. Solid lines represent stable fixed points - attractors (they correspond 
to the solution (4) with "— " sign before the square root); dashed lines represent 
unstable repellers (corresponding to "+" in (4)). 

We find two kinds of attractors: (i) unification (a = R when the leader's 
opinion wins, a = when it ceases to exist) and (ii) a stable cluster resulting 
from a solution of (4). In the SL-a space one can distinguish between three 
basins of attraction. Starting from a state in the area denoted as I, the time 
evolution leads to unification with a = 0. The stable cluster attractor divides 
its basin of attraction into the areas I la and lib. All states from III will evolve 
to unification with a = 20. Owing to the two possible signs of self-support 
parameter /3 in (4), the attractor and repeller are split. The space between 
their two parts enclose the "frozen" states that do not change in the course 
of time. These states correspond to local equilibria of the system dynamics 
similar to spin glass states. Thus, as a result of self-support, even repeller 
states can be stabilized. As one can see, the results of computer simulations 
fit the calculated curves very well. 

Taking into account the conditions (5), (6) and the two possible opinions of the 
leader one can draw a phase-diagram h — sl distinguishing the regions where 
different system states are possible [19,21]. Apparently, the system shows mul- 
tistability in a certain range of Sl and h. It depends on the history which of the 
states is realized, so we can observe a hysteresis phenomenon [19,21]. Moving 
in the parameter space Sl — Ji, while starting from different configurations one 
can have many possible scenarios of phase transitions [21]. 



2.3 Effects of social temperature 

It is obvious that the behaviour of an individual in a group depends not 
only on the influence of others. There are many more factors, both internal 
(personal) and external, that induce opinion formation and should be modeled 
somehow. In our model, we do this by means of a noisy dynamics, i.e. we use 
the equation (2) with the parameter T > 0. In the presence of noise, the 
marginal stability of unstable clusters due to the self-support is suppressed 
and they are no longer the stationary states of the system. The borders of 
the stable clusters become diluted, i.e. individuals of both opinions appear 
all over the group. Our simulations [19,21] prove that the presence of noise 
can induce the transition from the configuration with a cluster around the 
leader to the unification of opinions in the whole group. There is a well defined 
temperature T c that separates these two phases. To estimate the dependence of 
T c on other system parameters analytically, one can use a mean field approach, 
like methods developed in [19,21]. The two limiting cases of such an approach 
correspond to low temperature and high temperature approximations and are 
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discussed in the following. 



2.4 Low-temperature mean-field approximation 



For low temperatures T, i.e. for a small noise level, the cluster of leaders fol- 
lowers is only slightly diluted and its effective radius a(T) can be treated as an 
order parameter. One can then calculate the impact 1(d) acting on the group 
member inside (d < a) and outside (d > a) the cluster respectively [19]: 

I .(d) = -^-8aE(-^) +4^(4,^+2^-/3, (7) 



W) = S -j + 8a E (± , arcsin -\ - AR E (± , ^ + 2^ - P, (8) 



where E(k,ip) = Jq(1 — k 2 sin 2 a) l ^ 2 da is the elliptic integral of the second 
kind. Both functions are plotted in Fig. 2 for sl = 400. The system remains 
in equilibrium, therefore the impact on every individual is negative (nobody 
changes his/her opinion). It approaches zero at the border of the cluster which 
means that individuals located in the neighbourhood of that border are most 
sensitive to thermal fluctuations. We can however observe a significant asym- 
metry of the impact. It is considerably stronger inside the cluster. Individuals 
near the leader are deeper confirmed in their opinion, so they are also more 
resistant against noise in dynamics. When we increase the temperature start- 
ing from T ~ 0, random opinion changes begin. Primarily it concerns those 
near the border (the weakest impact). As a result individuals with adverse 
opinions appear both inside and outside the cluster. They are more numerous 
outside because of the weaker impact (cf. Fig. 2). 

Effectively, we observe the growth of a minority group. This causes the sup- 
portive impact outside cluster to become still weaker and the majority to 
become more sensitive to random changes. It is a kind of positive feedback. 
At certain value of temperature the process becomes an avalanche, and the 
former majority disappears. Thus, noise induces a jump from one attractor 
(cluster) to another (unification). Such a transition is possible for every non- 
zero temperature, but its probability remains negligible until the noise level 
exceeds a certain critical value that corresponds to our critical temperature 
T 

Using Eq. (2) and taking into account Eqs. (7) and (8) we can compute the 
probability Pr(<r = l)(r) that an individual at the distance r from the leader, 
shares opinion +1, which is assumed as the opinion of the leader. Then, the 
mean number of all individuals with opinion +1 may be calculated by inte- 
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grating this probability multiplied by the surface density (equaling 1) over the 
whole space: 



R 



n{a = 1)(T) = J Pi (a = l)(r)2nrdr. (9) 



This number equals the effective area of the circular cluster, so its radius is 



ln(a = l)(T) 



°co = v . ■ ( 10 ) 



7T 



The Eq. (10) is a rather involved transcendent equation for a{T) (it appears 
on the right hand side in ij(r) and / D ( r ))- For low temperatures T it has 
three solutions a(T) corresponding to a stable cluster, an unstable cluster and 
a social homogeneous state. The numerical solution for the radius of stable 
cluster is presented in Fig. 3 together with results of computer simulations. 
One should point out that the radius of the cluster a is an increasing func- 
tion of the temperature T for the reasons discussed above. At some critical 
temperature, a pair of solutions corresponding to the stable and the unstable 
cluster coincide [19,21]. Above this temperature, there exists only the solu- 
tion corresponding to the social homogeneous state. Fig. 4 shows the plot of 
the critical temperature T c obtained from (10) as the function of the leader 
strength sl together with results of computer simulations. 



2.5 High-temperature mean-field approximation 



For high temperatures or small values of the leader's strength sl, the cluster 
around the leader is very diluted and it is more appropriate to assume that 
there is a spatially homogeneous mixture of leaders followers and opponents, 
instead of a localized cluster with a radius a(T). It follows that at each site 
there is the same probability < p(T) < 1 to find an individual sharing the 
leaders opinion, and p{T) plays the role of order parameter. Neglecting the 
self-support {(5 = 0) one can write the social impact acting on a opponent of 
the leader at place x as [21]: 

I(x) = ^- + (2p-l)psJ D (x) + h (11) 



Jd( x ) = Sd r l/fi'dr — x|)rf 2 r is a function which depends only on the size of the 
group and the type of interactions. After a short algebra one gets the following 
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equation for the probability p(T) [21]: 



where I(x,p) is given by (11). Similar to equation (10) obtained for low tem- 
peratures, there are three solutions of Eq. (12): the smallest one corresponds 
to the stable cluster around the leader, the middle one to the unstable cluster 
which, in fact, is not observed, and the largest one to the unification. The size 
of the stable cluster grows with increasing temperature up to a critical value 
T c when it coincides with the unstable solution. At this temperature, a tran- 
sition from a stable cluster to unification occurs [21]. For T > T c , unification 
is the only solution, but it is no longer a perfect unification because due to 
the noise individuals of the opposite opinion appear. When the temperature 
increases further, p(T) tends to 1/2 which means that the dynamics is random 
and both opinions appear with equal probability. 



3 Modelling opinion dynamics by means of active Brownian par- 
ticles 

3.1 The model 

There are several basic disadvantages of the model considered in the previous 
chapter. In particular, it assumes, that the impact on an individual is updated 
with infinite velocity, and no memory effects are considered. Further, there is 
no migration of the individuals, and any "spatial" distribution of opinions 
refers to a "social" , but not to the physical space. 

An alternative approach [25] to the social impact model of collective opinion 
formation, which tries to include these features is based on active Brownian 
particles [22,?,?,?,?], which interact via a communication field. This scalar 
field considers the spatial distribution of the individual opinions, further, it 
has a certain life time, reflecting a collective memory effect and it can spread 
out in the community, modeling the transfer of information. 

The spatio-temporal change of the communication field is given by the follow- 
ing equation: 

d N 

— K(r, t) = J2 Si5 a , ai 8(r - r,) - 7 /i tr (r, t) + D h Ah a (r, t). (13) 
ot i=i 
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Every individual contributes permanently to the field h a (r,t) with its opinion 
<7j and with its personal strength Sj at its current spatial location Here, 5 a>ai 
is the Kronecker Delta, S(r — denotes Dirac's Delta function used for con- 
tinuous variables, N is the number of individuals. The information generated 
by the individuals has a certain average life time I/7, further it can spread 
throughout the system by a diffusion-like process, where Dh represents the dif- 
fusion constant for information exchange. If two different opinions are taken 
into account, the communication field should also consist of two components, 
a = {— 1, +1}, each representing one opinion. 

In this model, the scalar spatio-temporal communication field h a (r,t) [25], 
plays in part the role of social impact f used in [19,21]. Instead of a social im- 
pact, the communication field h a (r, t) influences the individual i as follows: At 
a certain location r,, the individual with opinion o~i = +1 is affected by two 
kinds of information: the information resulting from individuals who share 
his/her opinion, /i CT=+ i(i"j, t), and the information resulting from the oppo- 
nents /i CT= _i(rj, t). Dependent on the local information, the individual reacts 
in two ways: (i) it can change its opinion, (ii) it can migrate towards locations 
which provide a larger support of its current opinion. These opportunities are 
specified in the following. 

We assume that the probability pi(cri,t) to find the individual i with the 
opinion (Xj changes in the course of time due to the master equation (the 
dynamics is continuous in time): 

^Pi(a i: t) = ^w(a i \a' i )p i (a[,t) - p i {a i ,t)^w{a'^a i ). (14) 

where rates of transition probability are described in a similar way to Eq. (2) 
w(<7-|<7i) = ?7exp{[/v(rj,£) - K(ri,t)]/T} for a ^ a' (15) 

and iu(<Tj|<Tj) = 0. The movement of the individual located at space coordinate 
Yi is described by the following overdamped Langevin equation: 

+ yf2D^Si(t). (16) 



In the last term of Eq. (16) D n means the spatial diffusion coefficient of the in- 
dividuals. The random influences on the movement are modeled by a stochas- 
tic force with a 5-correlated time dependence, i.e. £(£) is white noise with 
(£i(t) £j(t')) = SijS(t — t'). The term h e (r,t) in Eq. (16) means an effective 
communication field which results from h a (r,t) as a certain function of both 



drj _ dh e (r,t) 
dt dr 
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components, h±i(r, t) [25]. Parameters on are individual response parameters. 
In the following, we will assume cej = a and h e = h a . 



3.2 Critical conditions for spatial opinion separation 



The spatio-temporal density of individuals with opinion a can be obtained as 
follows: 

N 

n a (r,t) = J Yj^wK* ~ r i )P(a 1 ,r 1 ...,a N ,r N ,t)dr 1 ...dr N (17) 



P(a, r, t) = P(ai, ri, o"at, rjv, t) is the canonical iV-particle distribution func- 
tion which gives the probability to find the N individuals with the opinions 
<7i, <7jv in the vicinity of ri, r/v on the surface A at time t. The evolution 
of P(gi,r,t) can be described by a master equation [25] which considers both 
Eqs. (15), (16). Neglecting higher order correlations, one obtains from Eq. (17) 
the following reaction-diffusion equation for n a (r, t) [23,?]: 



d 

n a (r, t) = - V [n a (r, t) a Vh a (r, t)] + D n An CT (r, t) 



dt 

Y [w(a'\a)n a (r,t) +w(a\a')n (T >(r,t) 



(18) 



with the transition rates given by eq. (15). Eq. 18 together with Eq. 13 form 
a set of four equations describing our system completely. 

Now, let us assume that the spatio-temporal communication field relaxes faster 
than the related distribution of individuals to a quasi-stationary equilibrium. 
The field h a (r,t) should still depend on time and space coordinates, but, 
due to the fast relaxation, there is a fixed relation to the spatio-temporal 
distribution of individuals. Further, we neglect the independent diffusion of 
information, assuming that the spreading of opinions is due to the migration 
of the individuals. Using h a (r, t) — 0, = s and D h = we get: 

h a (r,t) = -n a (r,t) (19) 

7 



Inserting Eq. (19) into Eq. (18) we reduce the set of coupled equations to two 
equations. 
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The homogeneous solution for n a (r, t) is given by the mean densities: 

n , N . . 

— where n = — (20) 



Under certain conditions however, the homogeneous state becomes unstable 
and a spatial separation of opinions occurs. In order to investigate these crit- 
ical conditions, we allow small fluctuations 5n a ~ exp (At + ikr) around the 
homogeneous state n a and perform linear stability analysis [25]. The resulting 
dispersion relations read: 



Ai(k) = -k 2 C + 2B ; X 2 (k) = ~k 2 C 

- (21) 



It follows that stability conditions of the homogeneous state, n a (r,t) = n/2, 
can be expressed as: 

T >T{= — , D>D c n = °™ (22) 

7 ^7 



If the above conditions are not fulfilled, the homogeneous state that corre- 
sponds to paramagnetic phase is unstable (i) against the formation of spatial 
"domains" where one of opinions a — ±1 locally dominates, or (ii) against the 
formation of a ferromagnetic state where the total numbers of people sharing 
both opinions are not equal. 

Case (i) can occur only for a systems whose linear dimensions are large enough, 
so that large-scale fluctuations with small wave numbers can destroy the ho- 
mogeneous state [25]. In case (ii), each subpopulation can exist either as a 
majority or as a minority within the community. Which of these two possible 
situations is realized, depends in a deterministic approach on the initial frac- 
tion of the subpopulation. Breaking the symmetry between the two opinions 
due to external influences (support) for one of the opinions would increase 
the region of initial conditions which lead to a majority status. Above a crit- 
ical value of such a support, the possibility of a minority status completely 
vanishes and the supported subpopulation will grow towards a majority, re- 
gardless of its initial population size, with no chance for the opposite opinion 
to be established [25]. 
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4 Conclusions 



This work discusses the possibilities of phase transitions in models of opinion 
formation which are based on the social impact theory (two opinions case). 
In the presence of a strong leader situated in the centre of a finite group, a 
transition can take place from a state with a cluster around the leader to a 
state of uniform opinion distribution where virtually all members of the group 
share the leaders's opinion. The transition occurs if a leader's strength exceeds 
a well defined critical value or if the noise level ( "social temperature" ) is high 
enough. The weaker the leader's strength is, the larger noise is needed. The 
value of the critical temperature can be calculated using mean field methods 
where either the existence of an effective value of the cluster radius (low tem- 
perature method) or a spatially homogeneous mixture of both opinions (high 
temperature method) is assumed. Numerical simulations confirm the analytic 
results. 

The extension of the social impact model is based on the concept of active 
Brownian particles which communicate via a scalar, multi-component commu- 
nication field. This allows us to take into account effects of spatial migration 
(drift and diffusion), a finite velocity of information exchange and memory 
effects. The reaction-diffusion equation for the density of individuals with 
a certain opinion is obtained. In this model, the transition can take place 
between the "paramagnetic" phase, where the probability to find any of op- 
posite opinions is 1/2 at each place (a high temperature and a high diffusion 
phase), the "ferromagnetic" phase with a global majority of one opinion (a low 
temperature and a low diffusion phase) and a phase with spatially separated 
"domains" with a local majority of one opinion (an intermediate phase). 
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Figure Captions 



Fig. 1. Cluster's radius a vs. leader's strength sl - phase diagram for circular social 
space. Interactions proportional to inverse of mutual distance (I oc 1/r). Lines 
correspond to analytical results, points to computer simulations. 

Fig. 2. Social impact I as a function of distance d to the leader. Leader's strength 
s L = 400. 

Fig. 3. Mean cluster radius a vs. temperature T; sl = 400. Results of calculation 
(solid) and computer simulation (dotted). 

Fig. 4. Critical temperature T c (above which no stable cluster exists) vs. leader's 
strength sl- Leader's opinion fixed (independent of the group). Line - calculations 
(Eq. (10)), points - simulations. 
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